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Nuclear gluon modifications are the least constrained component of current global fits to nuclear 
parton distributions, due to the inadequate constraining power of presently available experimental 
data from nuclear deep inelastic scattering and nuclear Drell-Yan lepton-pair production. A recent 
advance is the use of observables from relativistic nucleus-nucleus collisions to supplement the data 
pool for global fits. It is thus of interest to investigate the sensitivity of various experimental 
observables to different strengths of nuclear gluon modifications from large to small Bjorken x. In 
this work we utilize three recent global fits with different gluon strengths to investigate the sensitivity 
of three observables: nuclear modification factor, pseudorapidity asymmetry, and charge ratio. We 
observe that both nuclear modification factor and pseudorapidity asymmetry are quite sensitive to 
the strength of gluon modifications in a wide pseudorapidity interval. The sensitivity is greatly 
enhanced at LHC (Large Hadron Collider) energies relative to that at RHIC (Relativistic Heavy Ion 
Collider). The charge ratio is mildly sensitive only at large Bjorken x. Thus measurement of these 
observables in proton-lead collisions at the LHC affords the potential to further constrain gluon 
modifications in global fits. 

PACS numbers: 24.85.-|-p,25.30.Dh,25.75.-q 



I. INTRODUCTION 

Nuclear parton distribution functions (nPDFs) are an 
important ingredient in the theoretical description of 
nucleus-nucleus collisions using the framework of pertur- 
bative Quantum Chromodynamics (pQCD). Although 
the scale evolution of the nuclear parton distributions 
can be described perturbatively, nPDFs are intrinsically 
non-perturbative. This non-perturbative nature makes 
their theoretical determination from QCD hitherto in- 
tractable. Therefore, as obtains in the case of nucleon 
parton distribution functions (PDFs), the nPDfs are cur- 
rently most efficiently determined from global fits to ex- 
perimental data. Earlier global analyses relied heav- 
ily on fixed-target nuclear deep- inelastic scattering (DIS) 
and Drell-Yan (DY) lepton-pair production, with the at- 
tendant low precision and almost total lack of constraints 
on nuclear gluon distributions. The use of an extended 
data set, incorporating data on inclusive hadron produc- 
tion in deuteron-gold collisions, has been pioneered in 
with better constraints on gluon modifications. De- 
spite all these advances the nuclear gluon distribution is 
still the least constrained aspect of global fits to nPDFs. 

Thus questions related to the distributions of partons 
(in particular gluons) in nucleons and nuclei remain of 
current interest, both theoretically and experimentally. 
Of particular importance are the marked suppression of 
the nuclear modification factor at forward rapidities and 
the issue of gluon saturation (see, for example 0, Q). 
The Large Hadron Collider (LHC), even more than the 
Relativistic Heavy Ion Collider (RHIC), affords the op- 
portunity to study several important aspects of these nu- 
clear parton distributions in a center-of-mass energy and 



transverse-momentum domain where pQCD is expected 
to work well. For precision tests of pQCD applied to 
nucleus-nucleus collisions and also for various other ap- 
plications it is important to gain control over the uncer- 
tainties in the elucidation of nuclear gluon distributions. 
It is therefore beneficial to investigate the sensitivity of 
various observables to varying strengths of gluon modifi- 
cations in different Bjorken- a; regimes to further enhance 
the constraining ability of global fits. 

In general the description of nuclear collisions based 
on pQCD is a complicated task. Much of the compli- 
cation derives from the intrinsically complex nature of 
the nuclear environment in collisions of heavy nuclei. In 
addition to the initial-state modification of the PDFs, 
the cross section of high-px hadron production is influ- 
enced by final-state effects (such as jet energy loss) and a 
complicated geometry. To better understand the physics 
of pQCD in the nuclear environment, it is highly de- 
sirable to disentangle the different nuclear phenomena 
affecting high-pT {pt ^1-5 GeV) hadron production. 
In this respect, asymmetric light-on- heavy nuclear colli- 
sions, such as p(d)A collisions offer the distinct advan- 
tage of absence of final-state effects. In fact, one of the 
physics goals of the Run 8 at RHIC was to provide a 
high-statistics "cold" nuclear matter data set to estab- 
lish a definitive baseline for "hot" nuclear matter, cre- 
ated in gold-gold collisions. The benchmark role of the 
"the deuteron-gold control experiment" fo r e.g . energy- 
loss studies has often been emphasized [1, [l0| . Proton- 
nucleus and deuteron-nucleus reactions have also been 
used to study the Cronin-effect [ll|, [T^]. These asym- 
metric systems also offer unique information about the 
underlying dynamics not available in symmetric proton- 
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proton or nucleus-nucleus collision systems. For example, 
a new feature, offered by nonidentical colliding beams 
like d+Au, is the pseudorapidity asymmetry, examined 
in some detail by the STAR collaboration [IJl- In view 
of the importance of asymmetric collisions, we anticipate 
similar pPb program at the LHC. 

The strategy adopted in this study is as follows: we 
consider three global nuclear parton distribution sets 
0] with different strengths of the nuclear gluon modifi- 
cations. We are mainly interested in the sensitivity of 
observables to gluon modifications; therefore we consider 
these nPDF sets simply as sources of differing severity 
of gluon modifications. Further details about this aspect 
of the study are given in the next section. In order to 
avoid issues related to final-state effects which are present 
in nucleus-nucleus collisions, we limit our calculations to 
charged hadron production in both dcuteron-gold colli- 
sions at RHIC and proton- lead collisions at the LHC. 
Utilizing this framework, we investigate the sensitivity 
of three basic quantities directly observable experimen- 
tally: the nuclear modification factor, the pseudorapidity 
asymmetry, and the charge ratio, to different strengths of 
the gluon modifications. The calculated results are com- 
pared with available experimental data: nuclear modi- 
fication factor for charged hadrons from the BRAHMS 
Collaboration [l^ . and pseudorapidity asymmetry from 
the STAR Collaboration respectively. 

The paper is organized as follows: in Sec. |TT]we review 
the basic formalism of the coUinearly-factorized pQCD- 
improvcd parton model as applied to nucleus-nucleus col- 
lisions. This section also includes the definition of the 
minimum bias nuclear modification factor, pseudorapid- 
ity asymmetry, and charge ratio. We present the results 
of our calculations for these quantities in Sec. IIIII Our 
conclusion is contained in Sec. IIVI 



II. 



THEORETICAL FRAMEWORK AND 
EXPERIMENTAL OBSERVABLES 

A. Collinear factorization 



The minimum bias invariant cross section for the pro- 
duction of final hadron h from the collision of nucleus A 
(with mass number A) and nucleus B (with mass num- 
ber B) {A + B ^ h + X), can be written, in collinearly 
factorized pQCD, as 



d p 

2. d(j[ah^cd) 



'ibcd 



Dh/cizcQj) 



dt 



sS{s + i+u) (1) 



Here Xa and Xf, are parton momentum fractions in A 
and B, respectively, and is the fraction of the par- 
ton momentum carried by the final-state hadron h. The 



factorization and fragmentation scales are denoted as Q 
and Qf, respectively. As usual, s, t, and it refer to the 
partonic Mandelstam variables, and the massless parton 
approximation is used in this study as expressed in the 
delta function. The quantity da{ab^ cd)/di in Eq. ^ 
represents the perturbatively calculable partonic cross 
section, and Df^^i.{zc,Q'j:) stands for the fragmentation 
function of parton c to produce hadron h, evaluated at 
momentum fraction Zc and fragmentation scale Qf. For 
a recent discussion on the issue of collinear factorization 
in nucleus- nucleus collision see [l^ . 

The homogeneous nuclear parton distribution func- 
tions (nPDFs) FafA{x,Q^), for a nucleus with N neu- 
trons and Z ~ A — N protons, with A the mass number, 
can in general be written as 

Z N 

Fa/A{X,Q^) = ■jfa/p/A{x,Q'^) + -jfa/n/A{x,Q'^) 

= (1 - ^^)fa/p/A{x, Q^) + Ojfa/n/A{x, Q^) (2) 

with oj = N/A and < a; < 1. The function 
fa/p/A{x,Q'^) denotes the parton distribution for parton 
a in a bound proton in nucleus A and fa/n/Aix, Q^) the 
distribution in a bound neutron respectively, and are re- 
lated, for valence quarks, by isospin symmetry. They are 
expressible as convolutions of free nucleon parton distri- 
bution functions (PDFs) fa/N{x,Q'^) and a shadowing 
function Sa/Aix,Q'^) which encodes the nuclear modifi- 
cations of parton distributions in nucleus A. We em- 
ploy three different nuclear parton distributions in this 
study: the HKN07 [1] ("weak" gluon modifications), the 
EPS08 ("strong" gluon modifications), the EPS09 @ 
("moderate" gluon modifications) sets. Note that for 
HKN07 and EPS09 the qualifiers "weak" and "moderate" 
respectively are with respect to the central fits; inclusion 
of uncertainties generates considerable spread in the rel- 
ative measure of gluon modifications. For consistency 
reasons we use the underlying nucleon PDFs employed 
in the HKN07 nPDFs, namely MRS98 [H]. The resuffs 
are quite insensitive to the choice of nucleon PDFs in 
general. For the final hadron fragmentation we utilize 
the fragmentation functions in the DSS07 set The 
factorization scale as well as the fragmentation scale are 
set to evolve with the transverse momentum of the out- 
going hadron, according to Q — Qf — pT- Since pQCD 
calculations are generally not reliable at very low px , the 
lower limit on our calculations has been set to pT = 1.5 
GeV/c. With our scale choice, this corresponds to a min- 
imum Q2 (Q2-) of 2.25 GeV2, while EPS08, EPS09, and 

HKN07 give 1.69,1.69, and 1.0 GeV^ for minimum Q^, 
respectively. The minimum for the DSS fragmenta- 
tion functions is given as 1.0 GcV^. 

It should be noted that whereas both EPS09 and 
HKN07 are available in both leading order (LO) and 
next-to-leading order (NLO) parametrizations, EPS08 is 
available only in leading order. Our calculations are 
therefore to leading order. Since we consider charged 
hadrons, we need charge-separated fragmentation func- 
tions. These fragmentation functions arc also available 
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at LO and NLO, with the NLO set adjudged more ac- 
curate We have performed calculations using both, 
and the results presented here are for the NLO set. The 
implication of using NLO fragmentation functions in a 
LO calculation is treated in Jj| . 

A comprehensive treatment of the uncertainties on cal- 
culated results has not been undertaken in the present 
study. The overall uncertainties on the calculated re- 
sults obviously depend on three components: the nPDFs, 
truncation of the perturbative series for the hard par- 
tonic cross sections, and the charged hadron fragmenta- 
tion functions. AUhough both EPS09 and HKN07 have 
facilities for estimating the uncertainties in the nPDFs, 
we have not included these uncertainties in our calcula- 
tions. We have used the central fits in these two routines 
without recourse to their associated uncertainties. While 
the inclusion of uncertainties may be important for com- 
paring the relative efficacy of different nPDFs, the focus 
here is on the sensitivity of the observables to different 
gluon modification scenarios. It is thus more transparent 
for the purpose at hand to use the central fits in HKN07 
and EPS09, but bearing in mind the potential overlap of 
calculated results with the inclusion of uncertainties in 
these nPDFs. The calculated quantities in the present 
study are ratios, therefore truncation errors largely can- 
cel (for a comparison of LO and NLO nuclear modifica- 
tion factor, see, for instance, (H). The errors arising from 
fragmentation functions are basically unknown. 



(p(d) side) and backward rapidities (A side), it is of in- 
terest to study ratios of particle yields between a given 
pseudorapidity value (interval) and its negative in these 
collisions. Thus the pseudorapidity asymmetry, Y^^, is 
defined as the ratio of the yield in the backward ( "target- 
side" , negative rapidities) region relative to the yield in 
the forward ("projectile-side", positive rapidities) region: 

YUpt)=E,^ / . (4) 

The meaning of the pseudorapidity asymmetry is clear: 
a ratio greater than unity implies more yield in the back- 
ward (A side) region relative to the forward (p(d)) region 
while a ratio less than unity implies the converse. Since 
the nPDFs are an important ingredient in hadron pro- 
duction, the study of pseudorapidity asymmetry can offer 
valuable information in further constraining the nuclear 
parton distributions. As remarked above, the pseudora- 
pidity asymmetry is independent of a reference pp cross 
section and thus more "self-contained"; it is also more 
global in the sense that it encapsulates information on 
nPDFs at different Bjorken-a; values. For conciseness, we 
will refer to Y^'^ > 1 as positive asymmetry and Y^^ < 1 
as negative asymmetry. 

D. Charge ratio and rapidity asymmetry of charge 

ratio 



B. Nuclear modification factor 

The minimum bias nuclear modification factor com- 
pares, as a ratio, spectra of particles produced in nuclear 
collisions to a hypothetical scenario in which the nuclear 
collision is assumed to be a superposition of the appro- 
priate number of proton-proton collisions. The ratio can 
be defined as a function of pT for any produced hadron 
species h at any pseudorapidity rj: 



1 E,A''<j%/d''p 



(3) 



AB £;^d3a^p/d3p 

Nuclear effects manifest themselves in R\g{pT) values 
greater or smaller than unity, representing enhancement 
or suppression, respectively, relative to the pp reference. 

The determination of the nuclear modification factor 
involves the cross section for proton-proton collisions in 
the denominator, and is thus sensitive to the precision of 
the knowledge of the pp cross section. The other two ob- 
servables considered, the pseudorapidity asymmetry and 
charge ratio, are independent of a baseline pp cross sec- 
tion and therefore do not suffer from this constraint. 



Let us define the charge ratio as the ratio of the yield of 
negatively charged hadrons to the yield of the positively 
charged hadrons. Then we have 



Z\B{PT,ri) 
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E, 
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The charge ratio is expected to be more sensitive to the 
fragmentation functions than to the nuclear parton dis- 
tributions, especially as one moves toward forward (large 
positive) rapidities. Nevertheless, in the backward di- 
rection (negative rapaidites) and even near midrapidity, 
the charge ratio shows some sensitivity to nuclear par- 
ton distributions. This rapidity dependence is conve- 
niently displayed in terms of a pseudorapidity asymme- 
try of the charge ratio, Y^^{pt), defined as the ratio 
of the charge ratio at negative rapidities relative to that 
at positive rapidities. Technically this is equivalent to 
the ratio of the pseudorapidity asymmetry of negatively 
charged hadrons relative to the asymmetry of the posi- 
tively charged hadrons: 



vAB 



(pi 



Z\b{ptM 



Z\b{pt, 



(6) 



C. Pseudorapidity asymmetry 



III. RESULTS 



In general, as the mechanisms for hadron production 
in p(d)A collisions may be different at forward rapidities 



The gluon distribution is the least constrained, and 
therefore, constitute the largest source of differences be- 
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twcen the various global fits to nuclear parton distribu- 
tions (see [1, H^). Thus the differences in the predictions 
of the three sets under consideration stem largely from 
the differences associated with their respective giuon 
modifications. These differences are largest at the initial 
factorization scale {Ql = 1.69, 1.0 GeV^ for both EPS08 
and EPS09, and HKN07 respectively) and decreases with 
increasing due to evolution. With our choice of scales 
iQ = Qf = Pt) the results presented here are sensitive to 
these differences up to the highest pt considered, pr = 60 
GeV. 

Nuclear effects encoded in the nPDFs are a;-dependent, 
and thus different effects are present at different x: shad- 
owing {x < 0.1, depletion), antishadowing (0.1 < x < 
0.3, enhancement), EMC effect (0.3 < x < 0.8, deple- 
tion), and Fermi motion [x > 0.1, enhancement). Due to 
the Xfi-integration in Eq. ^ [x^miPT^'n) < Xh < 1), dif- 
ferent nuclear effects are superimposed, thus it is difficult 
to effect a puritanical isolation of these different nuclear 
effects. 

In the following subsections we present the results of 
our calculations for nuclear modification factor, pseudo- 
rapidity asymmetry, and charge ratio for the production 
of charged hadrons in dAu collisions at RHIC and pPb 
collisions at the LHC. It should be emphasized that the 
resuhs from both HKN07 and EPS09 nPDFs are from 
their respective central fits. 



Nuclear modification factor of charged hadrons 
at RHIC and the LHC 



1. Charged hadron modification factor at RHIC 

In Figure [T] we present our result for nuclear modi- 
fication factor of charged hadrons at RHIC. Data ex- 
ist for negatively charged hadrons at forward rapidities 
and for sum of charged hadrons around midrapidity, and 
the level of agreement of calculation with data is as 
shown in the figure. Our calculations extends over the 
positive and negative BRAHMS pscudorapidity intervals 
(—3.5 < rj < 3.5), and thus cover reasonably well the 
whole spectrum of nuclear effects, viz, shadowing, anti- 
shadowing, EMC effects and Fermi motion. As is well 
known, large backward (negative) rapidities correspond 
to large Bjorken-cc momentum fractions in the gold nu- 
cleus while large forward (positive) rapidities correspond 
to small Bjorken-.T momentum fractions. Thus as we 
move from backward to forward different nuclear effects 
come into play: from EMC/Fermi motion effects at large 
X to progressively stronger suppression due to nuclear 
shadowing at very small Bjorken-a;. A major limitation 
at RHIC is that the kinematics are such that the very 
small X region {x << 0.1) is not appreciably accessed, 
except at very small pt , where soft physics effects could 
be significant. This constraint is alleviated at LHC en- 
ergies where it is possible to access the low x even at 
relatively high px- 



A further complication is the presence of strong isospin 
effects on the nuclear modification factor of charged 
hadrons in deuteron-gold collisions, in particular at very 
forward rapidities. This is readily apparent from the 
appreciably different modification factors at forward ra- 
pidities for positively charged and negatively charged 
hadrons respectively, in contrast with the equality of 
these modification factors at forward rapidities in proton- 
nucleus collisions. These effects are somewhat minimized 
for the sum of charged hadrons, especially about midra- 
pidity (—1.2 < rj < 1.2), and so we therefore illustrate 
the influence of nuclear effects on the nuclear modifica- 
tion factor of sum of charged hadrons. 




5 10 5 10 5 10 5 10 5 10 5 10 5 10 

PylGeV] 



FIG. 1: (Color Online) Nuclear modification factor, RdAu, 
for charged hadron production in d-(-Au collisions at RHIC 
(v^sjviv = 200 GeV) from backward {rj = —3.5) to forward 
{ri = 3.5) pseudorapidities. The solid line represents the 
HKN07 nPDFs, the dashed line the EPS09 nPDFs, while the 
dot-dashed line represents the EPS08 nPDFs. The filled tri- 
angles denote the BRAHMS data [l||. 

Very backward rapidities (—3.5 < t] < —2.9): at the 
lowest Pt considered {pt ^ 1.5), a;^"„ ^ 0.2, and thus 
the modification factor is determined by the relative con- 
tributions from antishadowing, EMC effect, and Fermi 
motion, leading to RdAu > 1- As pT increases, 
becomes larger and therefore only the EMC effect and 
Fermi motion contribute, with RdAu < 1- Finally, at suf- 
ficiently high Pt, only the Fermi motion is present, thus 
leading to an enhancement {RdAu > 1)- 

Backward rapidities (—2.35 < 77 < —0.8): at very low 
Pt, a;^"„ is in the region of the advent of shadowing for 
—2.35 < 77 < —1.9; and already in the shadowing re- 
gion for —1.2 < 77 < —0.8 respectively. Thus there is 
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full contribution from antishadowing (enhancement) for 
—2.35 < 77 < —1.9 and some shadowing (suppression) for 

— 1.2 < f] < —0.8. At intermediate pxs there are con- 
tributions from both antishadowing and the EMC /Fermi 
motion. For the interval —2.35 < 77 < —1.9 the antishad- 
owing contribution becomes negligible for pr > 5 GeV, 
while it is still present though small even at high pT for 

— 1.2 < 7] < —0.8. At high pT the modification factor is 
influenced dominantly by the EMC effect. 

Midrapidities (—0.2 < 77 < 0.2): shadowing con- 
tributes appreciably at very low pT with negligible con- 
tributions at higher px- The physics at intermediate 
and higher px is determined by the relative contributions 
from antishadowing (particularly for pt < 7.5) and the 
EMC effect. 

Forward rapidities (0.8 < rj < 2.35): shadowing con- 
tributes significantly at very low pT for 0.8 < rj < 1.2 
and up till around pr ^ i GcV for 1.9 < rj < 2.35, al- 
though the effect is more pronounced at lower px- The 
modification factor at higher px is determined by a mix 
of antishadowing and the EMC/Fcrmi motion. 

Very forward rapidities (2.9 < rj < 3.5): here shadow- 
ing dominates at very low px and the effect persists up 
to ~ 4 GeV. The effect of antishadowing is seen in the 
rise of RdAu between 4 and 8 GeV while the EMC effect 
is responsible for the downward trend of RdAu at Pt > 8 
GeV. 

This simple analysis using the x-dependence of nuclear 
effects is of course modified for positively and negatively 
charged hadrons respectively. For negatively charged 
hadrons, the net effect is a rapid enhancement of the 
modification factor with px, especially at very backward 
and forward rapidities, and a less steep enhancement with 
Px around midrapidities. The converse is the case for 
positively charged hadrons, although the suppression is 
less steep with increasing px at both backward and for- 
ward rapidities. 

The observed predictions of the three different nPDFs 
sets reflect the different relative strengths of the nuclear 
effects encoded in these nPDFs, especially the nuclear 
gluon modifications. Even though there arc clear differ- 
ences between the predictions of EPS09/EPS08 relative 
to those of HKN07 across board, a reflection of strong 
versus weak gluon modifications, the distinction between 
EPS08 and EPS09 is not very apparent at the very back- 
ward and very forward rapidities, except at low px- In 
the interval —2.35 < rj < 2.35, due to different gluon 
antishadowing and EMC, the distinction is more appre- 
ciable, despite the similarity in shape of the modification 
factor. The kinematics at the LHC will permit going deep 
into the shadowing regime at forward rapidities, and thus 
one anticipates clear-cut distinctions in the predictions of 
EPS08 and EPS09 nPDFs for significant range of pr- 

We now compare our results with experimental data. 
The BRAHMS Collaboration [i| has measured nuclear 
modification factor for sum of charged hadrons in the 
intervals —0.2 < 77 < 0.2 and 0.8 < rj < 1.2, and for neg- 
atively charged hadrons in the intervals 1.9 < t] < 2.35 



and 2.9 < t] < 3.5. Figure [T] shows our result for both 
negatively charged hadrons and sum of charged hadrons 
for the relevant ry intervals. The effect of the different 
gluon modifications present in the three nPDFs is appar- 
ent: EPS08 consistently has the strongest suppression at 
very low px, while HKN07 has the weakest, for all pscu- 
dorapidity intervals. The effect of gluon antishadowing 
and EMC is also readily seen for both —0.2 < ry < 0.2 
and 0.8 < ry < 1.2. In the interval -0.2 < rj < 0.2, 
EPS08, due to its stronger gluon modifications, has the 
best agreement with data {px < 6 GeV) while HKN07 
has the least. All three nPDFs are in good agreement 
with data for both 0.8 < 77 < 1.2 and 1.9 < 77 < 2.35, 
although the reach of the data is quite limited for the 
latter interval. 

At very forward rapidities (2.9 < 77 < 3.5) the data 
show a suppression for the yield of negatively charged 
hadrons not adequately reproduced by calculations for 
Px > 2 GeV, albeit with the small reach in px- This 
suppression has generated a lot of interest, with diverse 
mechanisms being proffered to explain the phenomenon 
(sec, for instance, |20| - |23| . Since our purpose here is to 
compare predictions from different nPDFs, we will not 
go into discussions of the various explanations of the 
observed suppression. Some earlier studies of nuclear 
modification factor with different nPDFs can be found 

in [H!,!!!]!^. 



2. Charged hadrons modification factor at the LHC 

We now consider nuclear modification factor of charged 
hadrons in proton- lead collisions at the LHC. Figure [5] 
shows our results for positively charged, negatively 
charged, and sum of charged hadrons. Unlike what ob- 
tains in deuteron-gold collisions, the modification fac- 
tors for positively and negatively charged hadrons are 
the same, except at very backward rapidities (77 = —6 
and 77 = —4) where isospin effects in the Pb nucleus are 
appreciable, and the comments above applicable. The 
behavior of the modification factor across the spectrum 
of pseudorapidities mirrors the trend seen in d-f Au colli- 
sions. We illustrate this for the sum of charged hadrons 
at the very backward rapidity. 

At 77 = -6 and px = 1.5 GeV, xfJl^ ^ 0.07 and 
thus there is negligible contribution from shadowing. The 
modification factor thus reflects the competing contribu- 
tions from antishadowing/Fermi motion (enhancement) 
and the EMC effect (suppression). Higher transverse 
momenta translate to larger a;^^„ , with less contribution 
from antishadowing, leading to RpPb < 1- With px > 10 
GeV, the dominant contribution is from Fermi motion, 
with the attendant RpPb > 1. 

The major difference is the accessibility of the small- a; 
region even at backward rapidities and at relatively high 
transverse momenta for both midrapidities and forward 
rapidities. For instance, at 77 = —2 and px = 10 GcV, 
^min ^ 0.01, well within the shadowing region while at 
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r] = 2 and pr = 60 GeV, ~ 0.001, signifying a sub- 
stantial contribution from shadowing at this relatively 
high pt- In fact, at very forward rapidities ( ?/ = 4, 6) the 
influence of shadowing is so strong as to render Rppb < 1 
for all Pt considered. The strong sensitivity of the nu- 
clear modification factor to gluon modifications enables 
a clear distinction of the predictions of the three nPDFs 
sets utilized in this study. 



as a ratio of backward and forward nuclear modification 
factors [1^, the pp cross section being symmetric with 
respect to pseudorapidity. Thus it is rather straight- 
forward to apply the preceding analysis of the nuclear 
modification factor using the x dependence of nuclear ef- 
fects to the Pt dependence of pseudorapidity asymmetry. 
Pseudorapidity asymmetry has been studied for different 
asy mmetric systems and at different energies in p6l [28l - 
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FIG. 2: (Color Online) Nuclear modification factor, Rppt, for 
cliarged hadron production in p+Pb collisions at the L13C 
(\/sjv]v ~ 8.8 TeV) from bacltward (r; = —6) to forward 
{ri — 6) pseudorapidities. The solid line represents the HKN07 
nPDFs, the dashed line the EPS09 nPDFs, while the dot- 
dashed line represents the EPS08 nPDFs. The short dashed 
line is the result from a recent Color Glass Condensate (CGC) 
calculation [27| . 

The Color Glass Condensate (CGC) is an attrac- 
tive framework for addressing forward rapidity (small-x) 
physics (for a recent review see 0). For comparison pur- 
poses we show the predicted nuclear modification factor 
for charged hadrons at rapidities y = 2, 4, 6 from a recent 
calculation in this framework. Our result shows less 
suppression and a more dramatic pT dependence than 
the CGC. 



B. Pseudorapidity asymmetry of charged hadrons 
at RHIC and the LHC 

We now address pseudorapidity asymmetry in d-|-Au 
collisions at RHIC and p-|-Pb collisions at the LHC. A 
transparent way to understand the pt dependence of 
pseudorapidity asymmetry is to consider the asymmetry 



1. Asymmetry at RHIC 

In Figure [3] we present the pseudorapidity asymmetry 
for charged hadrons and sum of charged hadrons pro- 
duced in deuteron-gold collisions at RHIC. In the inter- 
val 0.8 < \ri\ < 1.2 the calculated asymmetry is practi- 
cally the same for both positively and negatively charged 
hadrons, and by implication, the sum of charged hadrons. 
For PT <6 GeV both EPS09 and EPS08 yield a positive 
asymmetry (> 1), with the EPS08 giving a larger asym- 
metry especially at low transverse momenta due to the 
stronger gluon modifications. Above pt — G GeV, the 
two nPDFs set predict a negative asymmetry (< 1) of es- 
sentially the same magnitude. The effect of a rather weak 
gluon modification in HKN07 is apparent: the predicted 
asymmetry at low pT is less than that of both EPS08 and 
EPS09. In fact, the predicted asymmetry is positive for 
all transverse momentum considered, with the asymme- 
try essentially unity at large pt- The result for the sum 
of charged hadrons is compared with the available exper- 
imental data from the STAR Collaboration [l^l for the 
interval 0.5 < \r]\ < 1.0. Since this interval is not far from 
midrapidity, the calculated asymmetry is rather weak, in 
agreement with experimental observation (3^ . [33j . 

Larger asymmetries are expected as we move away 
from midrapidity. This is readily apparent in the in- 
terval 1.9 < I77I < 2.35 where the calculated asym- 
metries for all three nPDFs sets are quite significant. 
While HKN07 predicts essentially the same asymmetry 
for both charge specie and the sum of charged hadrons, 
both EPS09 and EPS08 predict, at higher transverse 
momenta, a larger negative asymmetry for the posi- 
tively charged hadrons relative to the negatively charged 
hadrons. The effect of the different nuclear gluon mod- 
ifications in the three nPDFs is also apparent: EPS08 
yields the largest positive (negative) asymmetry at low 
(high) Pt while HKN07 yields the smallest, especially for 
positively charged hadrons and sum of charged hadrons. 
Also there are small but discernible differences between 
EPS09 and EPS08 for all transverse momenta, unlike the 
case for the interval 0.8 < I77I < 1.2 where both predict 
the same asymmetry for pT > Q GeV. Whereas both 
EPS09 and EPS08 yield positive asymmetry for pT < i 
GeV and negative for pT > ^ GcV, HKN07 gives negative 
asymmetry only for pT > 7 GeV. 

The above observations are also relevant for the asym- 
metries from the three sets in the interval 2.9 < \rj\ < 3.5. 
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2. Asymmetry at the LHC 
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FIG. 3: (Color Online) Pseudorapidity asymmetry, YdAu for 
charged hadrons at 0.8 < \r]\ < 1.2, 1.9 < |7?| < 2.35, and 
2.9 < \ri\ < 3.5 from d+Au collisions at RHIC (^/ilTiv = 200 
GeV). The solid line represents the HKN07 nPDFs, the 
dashed line the EPS09 nPDFs while the dot-dashed line cor- 
responds to the EPS08 nPDFs. The filled circles denote the 
STAR data [H for 0.5 < < 1.0. 



The difference between EPS09 and EPS08, though small, 
is quite distinct for all charged species and at all rele- 
vant pt- Both predict more significant asymmetry for 
positively charged hadrons than for negatively charged 
hadrons and for sum of charged hadrons. The physics 
of the asymmetry (a ratio) in this interval is governed 
largely by shadowing/antishadowing (denominator) and 
EMC/Fermi motion (numerator). Thus both EPS08 and 
EPS09 predict negative asymmetry for 2.8 < pr < 7.5 
GeV and positive everywhere else. On the other hand, 
HKN07 yields a positive asymmetry for px < 4.5 GeV 
and Pt < 8.5 GeV, with negative asymmetry for trans- 
verse momenta within these bounds. 

It is pertinent to remark at this stage that even though 
the predictions from EPS08/EPS09 on the one hand, and 
those of IIKN07 on the other hand are clearly distinguish- 
able, the kinematic reach at RHIC is such as to limit the 
clarity of the distinction between a strong gluon modifi- 
cation scenario as present in EPS08 and the somewhat 
more moderate gluon modification as in EPS09. The 
kinematics at the LHC allow such a distinction, and thus 
the study of pseudorapidity asymmetry at the LHC may 
shed some light on the actual strength of nuclear gluon 
modifications. 



Figure 2] depicts our result for pseudorapidity asymme- 
try of charged hadrons at the LHC. Let us consider the 
first interval: \r]\ =2. Even at this relatively small pseu- 
dorapidity the different predictions of the three nPDFs 
are clearly distinguishable. Both EPS09 and EPS08 pre- 
dict a negative asymmetry for px < 2 GeV, with the 
EPS08 yielding about twice the value of EPS09 at the 
lowest Pt considered. Thereafter both predict a posi- 
tive asymmetry which peaks around pT — 9 GeV and 
then decreases towards unity at very high pT- At this 
peak, EPS08 yields an asymmetry of about 30% while 
EPS09 yields about 20%. The asymmetry predicted by 
HKN07 is characteristically different: it rises approxi- 
mately monotonically as pt increases and is positive for 
all Pt considered. At the highest pT considered {pt = 60 
GeV), the asymmetry predicted by HKN07 is roughly 
the same as EPS08 for both negatively charged hadrons 
and sum of charged hadrons, and about 15% greater for 
positively charged hadrons. From Figure 0] it is apparent 
that the region 3 < pt < 20 GeV gives the best discrim- 
inatory ability for the considered pseudorapidity. 

The distinguishability of the different predictions of 
pseudorapidity asymmetry by the three nuclear parton 
distributions sets is appreciably magnified for the interval 
|?7| = 4, especially for pT < 20 GeV. In this pt region the 
predicted asymmetries are positive and approximately in- 
dependent of the charge specie, with EPS08 giving the 
largest asymmetry and HKN07 the smallest. Above this 
region the predicted asymmetries are positive for nega- 
tively charged hadrons and gradually turn negative and 
close in magnitude as px increases for positively charged 
hadrons and sum of charged hadrons. Thus the region 
Pt < 15 GeV affords the best opportunity to discriminate 
between the different scenarios of gluon modifications as 
enshrined in the three nPDFs. 

The above observations for jr;! = 4 apply equally well 
for the pseudorapidity asymmetry in the interval |r;| = 6. 
Here, due to phase space limitations, good discrimina- 
tory ability is limited to the transverse momentum re- 
gion pt < 10 GeV. While the predicted asymmetries re- 
main positive for all pT in the case of negatively charged 
hadrons, positively charged hadrons (sum of charged 
hadrons) exhibit a crossover from positive asymmetry to 
negative asymmetry around pT = 6 GeV {pt = 10 GeV) 
and a sharp rise back to positive asymmetry at around 
PT = 15 GeV. 

Thus contrary to what obtains at RHIC energies, pseu- 
dorapidity asymmetry at the LHC offers the potential of 
clearly distinguishing between various nPDFs with dif- 
ferent nuclear gluon modifications. 
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FIG. 4: (Color Online) Pseudorapidity asymmetry, Yppb for 
charged hadrons at = 2, \rj\ = 4, and \rj\ — 6 from p+Pb 
collisions at the LHC (-^sjvjv = 8.8 TeV). The solid line rep- 
resents the HKN07 nPDFs, the dashed line the EPS09 nPDFs 
while the dot-dashed line corresponds to the EPS08 nPDFs. 



Charge ratio and its asymmetry at RHIC and 
the LHC 



The top panel of Figure |S] shows the result of our cal- 
culation of the charge ratio in d+Au collisions at RHIC. 
The ratio is consistently below unity for all rj intervals 
and transverse momenta considered, with the implication 
of a suppression of negatively charged hadrons relative 
to the positively charged hadrons, in accord with expec- 
tations from fragmentation. The suppression is small- 
est around midrapidity where the charge ratio is close 
to unity, and becomes progressively more enhanced as 
one moves towards larger negative and positive rapidi- 
ties. In the very backward rapidity interval considered 
(—3.5 < T] < —2.9), the sensitivity of the charge ratio 
to nPDFs is most pronounced. While both EPS09 and 
EPS08 manifest practically the same behavior, IIKN07 
displays a greater suppression than both EPS09 and 
EPS08. This difference decreases significantly as we move 
towards positive rapidities, such that already at the in- 
terval 0.8 < ?7 < 1.2 all three nPDFs exhibit the same 
behavior for the charge ratio. Thus at very forward ra- 
pidities the charge ratio is independent of the differences 
in the nuclear parton distributions, and is thus solely de- 
termined by the fragmentation functions. 

The same trend is observed for p-|-Pb collisions at the 
LHC. Both EPS09 and EPS08 predict essentially the 
same behavior for all rapidities and transverse momenta. 



As in the case of d-|-Au collisions at RHIC, the sensitiv- 
ity is most pronounced for the very backward pseudo- 
rapidities and disappears at very forward rapidities. At 
midrapidity (?; = 0), all three nPDFs yield a ratio that 
is essentially unity. The HKN07 nPDFs set manifests 
quantitatively a greater suppression at very backward ra- 
pidities, same as exhibited in d-|-Au collisions at RHIC. 




FIG. 5: (Color Online) Top panel: Charge ratio, ZdAu, 
for charged hadron production in d-(-Au collisions at RHIC 
(v^sjvjv = 200 GeV) from backward (77 = —3.5) to forward 
{ri = 3.5) pseudorapidities. Bottom panel: Charge ratio, 
Zppb, for charged hadron production in p+Pb collisions at 
the LHC (^i]^ = 8.8 TeV) from backward (77 = -6) to 
forward (77 — 6) pseudorapidities. In general solid line rep- 
resents the HKN07 nPDFs, dashed line the EPS09 nPDFs, 
while dot-dashed line represents the EPS08 nPDFs. 

We now present the result for the pseudorapidity asym- 
metry of the charge ratio, . As remarked above, 
this is just the ratio of the usual pseudorapidity asymme- 
try of negative hadrons relative to the positively charged 
hadrons. Since the charge ratio is close to unity around 
midrapidity, the focus of interest is the behavior at both 
very forward (small (large) x in target (projectile)) and 
backward (large (small) x in target (projectile)) rapidi- 
ties. Thus the pseudorapidity asymmetry of the charge 
ratio helps to display the degree of asymmetry of the 
suppression relative to midrapidity. 

In Figure (5) we display the asymmetry for d+Au col- 
lisions (top panel) and p-|-Pb collisons (bottom panel). 
For deutcron-gold collisions, the asymmetry is positive 
and increasing with increasing 77 for both EPS09 and 
EPS08. In the case of HKN07, the asymmetry is con- 
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sistent with unity for the first two rj intervals, while it is 
slightly negative at large pr for 2.9 < Ir]] < 3.5. Let us 
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FIG. 6: (Color Online) Top panel: Pseudorapidity asymmetry 
of charge ratio, Vi*^" for charged hadrons at 0.8 < |7;| < 1.2, 
1.9 < |?7| < 2.35, and 2.9 < \r]\ < 3.5 from d+Au collisions 
at RHIC {^snn = 200 GeV). Bottom panel: Pseudorapidity 
asymmetry, Y^^*" for charged hadrons at|77|=2, |r;|=4, and 
|?7| = 6 from p+Pb collisions at the LHC {y^sJ^N = 8.8 TeV). 
In both cases the solid line represents the HKN07 nPDFs, 
the dashed line the EPS09 nPDFs, while the dot-dashed line 
represents the EPS08 nPDFs. 

now consider p+Pb collisions at the LHC. Both EPS09 
and EPS08 predict positive asymmetry for all rapidities 
considered, with the asymmetry increasing sharply with 
increasing rj. This trend is replicated by the HKN07 set. 
The degree of asymmetry is consistent for all three sets 
at 1 77 1 = 2, with slight differences at high pt for larger 
rapidities. 

IV. CONCLUSION 

In this study we have investigated the effects of differ- 
ent strengths of nuclear gluon modifications on charged 



hadron production at both RHIC and the LHC, using 
three recent nuclear parton distributions. We considered 
asymmetric light-on-heavy systems where final-state ef- 
fects (apart from fragmentations leading to final charged 
hadrons) are absent. The present treatment include ob- 
servables like nuclear modification factor, pseudorapid- 
ity asymmetry, and charge ratio for both negative and 
positive pscudorapiditcs and up to relatively high trans- 
verse momenta, thereby ensuring adequate coverage of 
all relevant nuclear effects encoded in the nuclear parton 
distributions. 

Each of the observables considered has its own 
strengths and weaknesses as a probe of nuclear effects in 
nuclear parton distributions. Nuclear modification fac- 
tor probes directly nuclear effects for any pT and at any 
pseudorapidity 77. Its drawback is the dependence on 
knowledge of the requisite pp yield. The associated mod- 
ification factor, central-to-peripheral ratio, while not de- 
pendent on the pp cross section, is geometry dependent. 
The pseudorapidity asymmetry, on the other hand, is in- 
dependent of the pp cross section, but requires the cross 
section of p(d) A collisions at both forward and backward 
rapidities. It is thus somewhat less "clean" than the nu- 
clear modification factor, but still has good potentials for 
constraining nPDFs. As shown above, the charge ratio is 
only sensitive in the backward region (large x in Au(Pb)); 
so it may offer some discriminating abilities for nuclear 
effects at large x. 

While there are some clear differences between the pre- 
dictions of the three nPDFs for the nuclear modification 
factor and pseudorapidity asymmetry at RHIC, the dis- 
tinctions are significantly enhanced at the LHC, even up 
to relatively high transverse momenta. In particular, at 
very forward rapidites where nuclear shadowing is ex- 
pected to dominate, the predictions of the nuclear modi- 
fication factor for all three nPDFs are directly correlated 
with the strengths of their gluon modifications. This is 
also true for the predicted pseudorapidity asymmetries. 

In conclusion, asymmetric p(d)A collisions, especially 
at the LHC, have great potential not only to shed light 
on the underlying dynamics of rclativistic nuclear colli- 
sions, but also as a source of data to help better constrain 
nuclear gluon distributions. 
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